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We investigate collective flavor oscillations of supernova neutrinos at late stages of the explosion. 
We first show that the frequently used single-angle (averaged coupling) approximation predicts 
oscillations close to, or perhaps even inside, the neutrinosphere, potentially invalidating the basic 
neutrino transport paradigm. Fortunately, we also find that the single-angle approximation breaks 
down in this regime; in the full multiangle calculation, the oscillations start safely outside the 
transport region. The new suppression effect is traced to the interplay between the dispersion in 
the neutrino-neutrino interactions and the vacuum oscillation term. 

PACS numbers: 97.60.Bw, 14.60.Pq 



1. Introduction. - It was first suggested more than 20 
years ago [J [5] that interactions between neutrinos in a 
core-collapse supernova could impact their flavor evolu- 
tion. Yet, modeling this effect has proven very challeng- 
ing. Although recent years have brought considerable 
progress, we still do not have a complete picture of this 
phenomenon and surprises continue to emerge. 

In its most general form, the problem involves model- 
ing of a coupled neutrino ensemble, with refraction and 
scattering, in an anisotropic and fluctuating environment 
of an exploding supernova (e.g., [3])- This is a formidable 
task, even for today's fastest supercomputers. Yet, it is 
also unnecessary if flavor transformations always happen 
safely outside of the neutrinosphere. In that case, the 
full quantum-kinetic equations could be replaced by the 
flavor oscillation equations for streaming neutrinos [H [S] , 
which only keep the coherent (forward scattering) part of 
the full Hamiltonian [BHTT] . Modern analyses begin with 
this assumption and check its consistency a posteriori. 

Assuming additionally that the effect can be modeled 
in spherical symmetry, the problem becomes manageable, 
even though still computationally intensive: with O(IO^) 
energy bins and 0(10^) angular bins one needs to solve 
millions of coupled differential equations. This task was 
first tackled only in 2005, on a supercomputer [UlIIS], 
revealing rich, unexpected physics: neutrinos leaving the 
"self-coupling" region undergo large flavor transforma- 
tions and settle into characteristic nonthermal "split" 
spectra. Crucially, the transformations in [T^HIl] in- 
deed occurred safely outside the neutrinosphere (starting 
at 70-80 km), providing a highly non-trivial consistency 
check for the free-streaming framework. 

Are further simplifications of this problem possible? In 
this letter, we will focus on what is known as the single- 
angle approximation [f51 116j . In this approximation, one 
assumes that neutrino wavefunctions are independent of 
neutrino trajectories and that it is sufficient to follow 
the neutrino flavor evolution along a single, representing 
trajectory. The full calculation, as described above, is 
referred to as multiangle. 

In [12], the single-angle and multiangle numerical re- 
sults were observed to be quite similar. The same obser- 
vation was also later made in [T3] , where it was suggested 



as a general result. This came as a very welcome devel- 
opment for the field, as single-angle calculations allow 
one to quickly explore various regimes of collective oscil- 
lations. In recent years, this framework has been used to 
uncover more rich physics, for example, the interplay of 
collective oscillations with the usual MSW effect [TTl [TB] , 
the possibility of multiple spectral splits [T^] , the appear- 
ance of the "mixed" spectra and three-flavor instabilities 
PUI \7T\ and other effects (see, e.g., [22] for a review). 
With so much staked on the validity of the single-angle 
approximation, it is of great interest to establish any lim- 
itations it might have, and physical reasons for them. 

2. Setup of the calculations. - For reasons that will be- 
come clear shortly, we choose the conditions that occur 
several seconds into the explosion. We adopt the late- 
time spectra of Ref . [IS] , as previously used in [TH HO] 
[36) and take the density profile of the neutrino-driven 
wind to be p = po{W km/r)^, with po = 10^ g/cm~'^, 
and Ye = 0.4, inspired by [24]. For definiteness, we start 
by choosing the inverted hierarchy (IH) of the neutrino 
masses, with 613 — 0.01 and the standard oscillation pa- 
rameters, as in [2Dj. The neutrinosphere is modeled as 
a sharp sphere with i? = 10 km, without limb darken- 
ing |12j . These ingredients represent a rather standard 
reference setup in the field. Assuming the consistency of 
the framework, the neutrinos are held unoscillated "by 
hand" for radial distances less than rgtart ~ 40 km (see 
later). We solve the multiangle equations (here and in 
what follows) using the same numerical code FLAT [5S] 
that was used for 3-flavor multiangle calculations in |26j . 

3. Single- vs. multi-angle results. - The resulting evo- 
lution is shown in Fig. [Tj Let us begin by considering the 
single-angle curve (marked "single-angle"). The quan- 
tity plotted captures the deviation of the spectrum from 
the initial one: |(<^-P)| ^ 0(1) indicates large transforma- 
tions, while |((5-P)| ^ 1 means the spectrum is nearly un- 
changed [37] . A striking aspect of this result is that the 
neutrinos are unstable to oscillations the moment they 
are released. What this suggests is that the oscillations 
could be present close to the neutrinosphere. Indeed, re- 
peating the calculations with smaller r-gtart , such as 30 km 
or even 20 km, we again find that the system is immedi- 
ately unstable to oscillations. This finding is independent 
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FIG. 1: (Color online) The onset of the collective oscillations 
in the multiangle and single- angle cases. The onset is insen- 
sitive to the background density, but depends strongly on the 
neutrino luminosities (as marked). 
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FIG. 2: (Color online) Same as Fig.[T] but with no matter and 
only two flavors. The onset of the oscillations and its depen- 
dence on the neutrino luminosities are virtually unchanged. 



of the exact single-angle prescription: we found immedi- 
ate oscillations whether we used the coupling for the ra- 
dial trajectory, for the tangential, or for the one emitted 
at 45° at the neutrinosphere. Thus, the short-distance 
suppression observed in the single-angle calculations in 
[T2HI1] is not a general result, but a consequence of the 
specific choice of conditions (emitted spectra) [55] , 

Taken at face value, this would mean neutrino trans- 
port has to be modeled together with collective oscilla- 
tions. Fortunately, the multiangle calculations come to 
the rescue. They are also shown in Fig. [T] by the thick 
curve. In this case, the oscillations are suppressed close 
to the protoneutron star. The flavor instability is seen to 
develop only at ri^st ^120 km. What we have, therefore, 
is a failure of the single-angle approximation and not of 
the neutrino transport paradigm. 

4. Deconstructing the effect. - Let us investigate the 
origin of this suppression. There are two conditions un- 
der which the single-angle approximation has been found 
to fail: (i) when the number fluxes of and Pe are highly 
symmetric [27l [28] , or (ii) in very dense matter [29l [30] . 
Neither of these conditions applies here. In particular, 
the high-density criterion is rig- — 71^,+ > |29| . This 
condition could be present during the initial stage of the 
explosion, when dense matter is piled up above the neu- 
trinosphere. It should disappear, however, once the shock 
is pushed out. (This is the reason why we here focus on 
the late-time conditions.) 

We can see explicitly that matter does not play a role 
if we increase and decrease the background density by a 
factor of five: as the dashed and dash-dotted curves in 
Fig. [T] illustrate, the instability still sets in at ^ 120 km. 

In fact, we can remove background matter altogether. 
In Fig. |2] we show the results of the calculations for such 
a simplified setup. The thick curve shows that the large 
conversion still occurs around ri^st ~ 120 km. The cal- 
culation in this case is done with only two vacuum eigen- 



states separated by the atmospheric splitting, Am^^^. 
Thus, the three-flavor effects are not playing a crucial 
role in the location of the onset of the instability. 

The only remaining ingredient in this simplified sys- 
tem is the V — V interactions. These are different in the 
single- and multi-angle calculations: in the latter case, 
the strength of the potential felt by a neutrino depends on 
its direction. This dispersion of the interaction strength 
must be the cause of the suppression. 

5. Discussion. - First, let us review how dense mat- 
ter suppresses collective modes. In spherical symmetry, 
collective oscillations develop radially, so that on any ra- 
dial shell the oscillation phase is the same for all neutri- 
nos. Yet, in the presence of matter and without v — v 
coupling, neutrinos on different trajectories would accu- 
mulate different phases between radial shells. Indeed, 
a neutrino making an angle with the radial direction 
travels 8r j cos •& between two spherical shells separated 
by br. (Unlike the common matter potential, this dis- 
persion, AHj^s^t ^ ^^mat(l/cosi?max " 1), cauuot bc re- 
moved by using the corotating frame technique in j31j.) 
The V — V interactions need to overcome Ai/n^at to keep 
the neutrinos locked in a collective mode, resulting in the 
criterion n^- — ng+ < n^, [29l [30] for the oscillations. 

It is easy to quantify how the potential felt by a neu- 
trino depends on the angle {) its momentum makes with 
the radial direction. The interaction strength between a 
pair of neutrinos making an angle is proportional to 
(1 — COS0). Integrating over the bundle of intersecting 
neutrino rays, one easily finds (c/. [TO] ) 

H,, ^ V2GF[n,Ar) ~ n^A^W ~ cosz9(l + x)/2], (1) 

where x = yj\ — B? jr"^ and cosi9 — y^l — sin^ dRB? jr"^ 
in terms of the angle of emission "Qb.- In turn, the number 
densities are related to the fluxes /j = Li/ {Ei) by 
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FIG. 3: (Color online) (Top): The pattern of flavor trans- 
formations at the start of the oscillations, corresponding to 
the simulation shown by the thick curve in Fig.[l] (Bottom): 
Isocontours of H^^ + if^ac - {H'(^r + HZ'')/2, in units of 
Am^tm/(20 MeV). Both are plotted for r = 133.6 km. 



Strictly speaking, one should also include the geomet- 
ric factor l/cosi?, just like dense matter. Its effect, how- 
ever, appears only at a higher order: for sufficiently large 
r, 77^1/ oc /Ar^ -f 19^/2 -f while the geometric cor- 
rection appears only at 0(i9^). This illustrates an impor- 
tant point: the dispersion in AiJmat is small compared 
to i?niat , while the dispersion in Hyi, is of the same order 
as H^ij. Therefore, the latter effect is more important, 
unless the matter is very dense. 

To understand what sets the location of the instabil- 
ity J'inst, let us first recall that in a single-angle calcu- 
lation collective transformations are driven by the vac- 
uum term -ffvac — lS.m? /2E^, or more accurately, by its 
dispersion AH^i^c with energy (a common part can be 
removed, just like for the matter potential [21]). In the 
multiangle setup, the dispersion of H^i, with angle t?, 
AHiJi, ^ Hvi,(d = arccosx) — H^^^d = 0) adds in. Since 
A.Hi,y (just like Hyy) scales as r~^, it dominates over 
AiJvac at small r. The regime when AH^^ "overpowers" 
AHvac is seen to give no oscillation. The flavor conver- 
sion therefore starts when the two become comparable. 



AH, 



(3) 



Plugging in the numbers, one indeed finds rjnst ^ 120 km. 
Furthermore, since AHi,, falls off as r^'', while ifvac is 
independent of r, ri^st should depend on the luminosity 
as lI,^^. This is precisely what is seen in Figs, [ijandji] 

Further insight into this physics can be obtained from 
analyzing the oscillation mode in the {Ei,^^ji) plane. In 
Fig. |3j in the top row we show the pattern of conversion 
at r — 133.6 km, where the instability just triggered con- 
version. As we have argued, the evolution at this point 
should be driven by a combination of AiJvac and AH^y. 
In the bottom rows of Fig. [3] we show the isocontours of 
We can explicitly see that the two sources 



of the dispersion are comparable at this radius [32] and 
furthermore, that the pattern of conversion matches the 
pattern of the isocontours remarkably well. 

We also investigated the normal mass hierarchy (NH) 
scenario. The oscillations in this case for our reference 
model also start at r-jnst ^ 120 km. We confirmed that 
Hnst again scales with neutrino luminosity as L^/^. The 
oscillation pattern in the {E^^dp) plane is different from 
the IH case. Yet, once again, it follows the isocontours of 
Hvac + remarkably well. (In this case, the left and 
right panels of the bottom row in Fig. [3] switch places, 
because Ann? has the opposite sign.) 

Animations depicting the fiavor evolution as a function 
of r, for both IH and NH, are available online [40] . 

6. Final neutrino spectra. - In Fig. [4] we plot the 
spectra "at infinity" for both neutrino mass hierarchies. 
Compared to the single-angle calculations, the multian- 
gle ones give features that are somewhat smeared out, 
especially in the NH case, but are qualitatively similar 
|41) . This behavior, however, is not general. Our investi- 
gations show that for other spectra, the single- and mul- 
tiangle calculations give qualitatively different answers 
not just at intermediate radii, but also at infinity. These 
results will be reported elsewhere. 

7. Generalizations. - For applications, one needs to 
know how the suppression effect manifests itself in a vari- 
ety of conditions that may exist in a supernova. We have 
already seen that ri„st is insensitive to the matter pro- 
file, but varies with the overall neutrino luminosity. An- 
other obvious factor is the radius of the neutrinosphere 
R, which decreases with time. Since AH^^ scales as 



cx R , while also being proportional to r 
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should be proportional to v i? 
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FIG. 4: (Color online) Spectra of (left) and (right) at in- 
finity, for both IH (top) and NH (bottom). Multiangle (M.A.) 
results are shown with thick solid curves, single-angle (S.A.) 
with filled regions. Initial spectra are also shown, as marked. 
Turbulence (32, and shock front [33] effects are ignored. 
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We have verified that doubling R shifted rinst from 120 
km to 170 km. 

Less trivial is the dependence of the effect on the neu- 
trino spectra. We have observed that, as the spectra are 
varied, the coefficient implicit in Eq. ([3| and the starting 
oscillation pattern in the {E^^-dji) plane change. As al- 
ready mentioned, for sufficiently different spectra (small 
Vx fluxes), oscillations at small r could be suppressed by 
the known single-angle mechanism [31) . This was the 
regime studied the original multiangle calculations [12] 
and in many subsequent papers. The study of different 
initial spectra is forthcoming. 

8. Conclusions. - We have described a new effect 
that suppresses collective neutrino oscillations close to 
the protoneutron star. The effect was traced to the dis- 
persion of the neutrino-neutrino interaction strength. A 
physical condition behind this suppression was discussed. 

The suppression has implications for supernova mod- 



eling: with collective flavor oscillations suppressed at the 
neutrinosphere, it is not inconsistent to treat the trans- 
port of neutrinos separately from collective oscillations. 
The suppression also impacts the nucleosynthesis yields, 
particularly the r-process, where the single-angle approx- 
imation gives a grossly inaccurate answer [34 . 

The physics of collective flavor transformations contin- 
ues to surprise us with its richness. As more conditions 
are investigated, a more complete physical understand- 
ing of the phenomenon will emerge. The role of numerical 
calculations in this field, as a tool of scientific discovery, 
has been indispensable. 
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work was supported by the DOE Office of Science and 
the LANL LDRD program. 
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